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STRESS ANALYSIS BY RECURRENCE FORMULA OF REINFORCED 
CIRCULAR CYLINDERS UNDER LATERAL LOADS 
By John E. Dub erg and Joseph Kempner 

SUMMARY 


A recurrence formula is developed for the Btress analysis of 
reinforced circular cylinders loaded in tile planes of their rings • 

In contrast to the elementary engineering analysis, deformations 
of rings and sheet are considered. The recurrence formula in 
conjunction with appropriate boundary equations can be used to 
obtain sets of simultaneous linear algebraic equations. The solu- 
tions of these equations enable the stress analyst to find the 
shear flows and direct stresses in the shoot, as well as the 
shear forces, axial forces, and bending mcments in the rings. 

In order to reduce the amount of computation involved in the 
stress analysis of relatively long reinforced cylinders, an approxi- 
mate method of analysis is presented. In this method the cylinder 
under consideration is assumed to be infinitely long, and the 
recurrence formula Is treated as a fourth-order finite-difference 
equation. It is recommended that the approximate solution be 
utilized for the stress analysis of cylinders loaded at rings 
located two or more bays from external restraints. 


INTRODUCTION 


Experimental data on stresses in reinforced circular cylinders 
indicate the inadequacy of the elementary theory of bending and 
torsion when applied to the relatively flexible shell structures 
used in airframe construction. Several investigators have pre- 
sented methods for the stress analysis of cylinders laterally loaded 
at the reinforcing rings (references 1 to 3)* The theory of refer- 
ence 1, developed only for the case of a one bay cylinder, involves 
the assumption that stringer strains can be entirely neglec'ted and, 
consequently, leads to inaccurate results. The more precise theory 
of references 2 and 3, developed for cantilevered cylinders having 
identical bays, becomes tedious and unwieldy when extended to non- 
uniform cylinders. 
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The present paper contains the development of a general 
recurrence formula suitable for the stress analysis of cylinders 
that may be nonuniform in construction,, arbitrarily supported at 
the boundaries, and arbitrarily loaded in the planes of the rein- 
forcing rings . The development is based upon tho maintenance of 
continuity of deformation between the rings and shell . In any 
particular problem the recurrence formula together •with appro- 
priate boundary equations are used to obtain sets of simultaneous 
linear equations for the corrections to tho stresses given by the 
elementary theory. (For a cantilevered uniform cylinder tho results 
obtained in this manner are identical with those obtained by tho 
method of reference 2 or 3 . } 

If a cylinder is composed of many bays, as in conventional 
fusolago construction, the number of simultaneous equations 
requiring consideration may bo .prohibitive • For a uniform cylinder, 
however, good approximations to tho correction strossos can be 
obtained if tho cylinder is assumed to be infinitely long. The 
recurrence formula for this case is solved as a homogeneous finite dif- 
ference equation of the fourth order and yields a relatively simple solu- 
tion. For practical purposes this solution can be applied to 
arbitrarily supported cylinders provided the loads are located a 
few bays from external restraints. 'When tho recurrence formula, 
together with the boundary equations presented, is applied to a 
cantilevered uniform cylinder discussed in reference 3, good agree- 
ment is obtained among the recurrence -formula solution, difference- 
equation solution, and experimental stresses. 


SYMBOLS 


IL 3 

B = Et ' r2 
G-tL 2 

C function of ring loading 



1 


Young's modulus 
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G 

H 

I 

L 

M 

M c 

P 

Q 

B 

T 

V 

a, b 
c 

1, k 

m 

n 

4 

t 

t* 

U, V, w 
7/ Z 


shear modulus 

axial force in ring 

moment of inertia of cross section 

length of bay 

bending moment 

concentrated ring bending moment 
radial load 

statio moment about neutral axis of cross-sectional area 
lying between extreme fiber and plane under considera 
tion 

radius of cylinder and ring 
tangential load on ring 
Bhear force 

Fourier coefficients in Fourier expansions of q. 

distance from neutral axis 

general numbers of bay or ring 

designation of root bay ! ' . 

general number of Fourier coefficient 

shear flow in skin 

thickness of skin 

effective sheet thickness, that is, thickness of all 
material carrying bending stresses in cylinder if 
■uniformly distributed around perimeter 

axial, tangential, and radial displacements of points 
on cylinder 

axial, tangential, and radial coordinates of cylindor 
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a 


arbitrary constant of integration 


i. - 3 * sLls 

n 3A 7 


n 


2 - 6B 


7 = -2 + 

n 12A7 



4 j (i, 15 . constants dependent upon bay lengths 



a longitudinal direct stress in Bkin 

cp angular coordinate of point on cylinder 



Subscripts: 

B rigid 


m moment . 

r radial 

t tangential 
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STRESS ANALYSIS OF REINFORCED CYLINDERS 
Inadequacy of Elementary Theory 


The elementary engineering theory for "bending and torsion of 
reinforced cylinders loaded at the ring reinforcements yields the 
well-known formulas Mc/l for direct "bending stress, VQ/lt for 
shear stress due to "bending, and T/2At for shear stress due to 
torsion (where T and A are the torq.ue on cross section and the 
area inclosed "by perimeter of cross section, respectively) • This 
simple theory is based upon the assumption that radial displace- 
ments of both rings and sheet can be neglected. Since the dimen- 
sions of most mcnocoque structures are ouch that radial displace- 
ments of the structural components cannot be ignored without 
appreciable inaccuracies the results of analysis, the elementary 
theory must be modified 60 es not only to satisfy the laws of 
statics but also to maintain continuity between rings and sheet. 

The present development, consequently, is directed towards finding 
self -equilibrating stress distributions that, when superimposed 
upon the elementary stress distributions, yield results which, in 
addition to satisfying the lavs of statics, preserve the continuity 
of the structure. Tbese correction stresses are found from the 
recurrence formula that is developed herein. 


Basic Assumptions of Present Theory 

In the development of the recurrence formula that can be 
used to obtain the desired stress corrections, several simplifying 
assumptions are made. That part of the sheet area which is con- 
sidered to resist normal stresses is added to the stringer area 
and the combination is uniformly distributed about the periphery 
of the cylinder. This resulting combination is an effective sheet 
thickness t-' that resists normal stresses. The actual Bheet area 
is considered capable of supporting only shear stresses . It then 
follows that within a bay the shear stresses vary in the circum- 
ferential direction but are constant in the longitudinal direction. 
Inexten&ional deformation of rings and sheet is also assumed, and 
Poisson's ratio is considered to be zero. 


Development of Recurrence Formula 

Procedur e.- For the skin of any bay i of the structure (see 
figs . 1 and 2) , the corrections to the elementary shear flow, 
direct stress, axial displacement, and radial displacement are 
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each expressed as Fourier series with undetermined Fourier coef- 
ficients. Through static, elastic, and - geometric considerations 
of rings and sheet, a recurrence formula is obtained relating the 
Fourier coefficient of the shear flow of any bay i with the coef- 
ficients of the two bays on each side of bay i, that is, bays i+1 
and i+2 and bays i-1 and i-2. From the recurrence formula, simul- 
taneous equations may be obtained from which the values of the 
shear-flow coefficients are determined. With these values the 
loads and stresses in the rings and sheet can be found. 

Sheet s tresses and deforma tions.- The system of coordinate 
axes to be used is shown in figures 1 and 2. Positive displace- 
ments in x-, y-, and z-directions are designated u, v, and v, 
respectively. For convenience, the external loading on the rein- 
forcing rings of a cylinder is considered to be either symmetrical 
or autifijrametrical about cp = 0°. (Soa figs. 1 and 2.) 3h accord- 
ance with the basic assumptions the corrections to the elementary 
shear flow, direct stress, axial displacement, and radial displace- 
ment at any point cpj in bay i can be expressed for symmet- 

rical loading as the Fouriop expansions 


^ e ln sin nq> 
xt=2 


(la) 


CO 

a i( x i' *p) " *in( x i) cos «p 


(lb) 


rt=2 


CO 

V 1 


\ (x ± , 9) = l_ cos «p 


(1c) 


rt=2 


co 

\( X V *) - 2_ V m( 3 i) 003 “P 

rt=2 


(Id) 
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respectively, in which a in , tr in( 3 ^.)> u in( x i)> and v in(*i)- "*» 

Fourier coefficients. Inasmuch as only corrections to the ele- 
mentary stresses and displacements are desired, Fourier terms 
corresponding to n = 0 and n = 1 are omitted since they corre- 
spond to the elementary stress and. displacement distributions • ; 

If antisymmetrical loading is considered, the harmonic func- 
tions in equations (1) are replaced -hy- their ' c'Ofuhctions. It is 
then convenient to designate the Fourier coefficient of the shear 
flow "by tj_ n . 

Selationship s among sheet stresses and deformations . - Within 
any hay i the following relationships exist, (fig. 2): hy the- . 

equilihrium equation ' 

: . : . v , . . t , a-afe,. <p) + i Mil o : - ' ( 2 a) 


hy Hooke's law for direct -stress 


.(*!, CP ) 5=: 


E 


5u i (x 1 , . t o) . ■ ... 


‘(2b) 


.hy Hooke's law for' shear stress 


9.1(9) 1 ° U j 9) 5? j(-i^. ?) . , .. 

■ G-t^ B ’* Sep "bx i 


(2c) 


f. . . \ 


and hy the inexbensional deformation equation (p* 208 of 'refer- 
ence ii-) 


— *) - 0 


(2d) 
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•where 

t^ actual skin thickness of bay i 

t effective skin thickness of bay i 

B radius of cylinder 

E Young's modulus 

G shear modulus 

v^^, circumferential displacement of any point in bay i 

If equations (la) and (lb) are substituted into equation (2a) 
and if coefficients of like cosine terms are equated, the following 
expression for the Fourier coefficient a in (xj\ is obtained: 



Integration of this equation yields 


. nxj_ 

or in( 3! i) = " £7T a in + ff in(°> (3) 

in Ttdilch d^ n (0) is the direct-stress Fourier coefficient at x^mO. 

Similarly, elimination of c ir ^x i ^( from equations (lb), (1c), 
(2b), and ( 3 ) and subsequent integration gives 

M 11 ) = ' + r , ^ <0> + “ to<0) (4) 

is the axial displacement coefficient at x^ = 0. 


in which u^fO) 
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Elimination of v^x^, <p^ from equations (2c) and (2d) yields 

5yj(xi ? cp) i 5q 1 (y) i 

oxj Gt i Sep E ^2 

Substitution of equations (la), (1c), (id), and (4) into this 
relationship and integration yields the following expression for 
the radial displacement coefficient: 


v. 


in 


( z i) - 


racj , 

Gt, 


in 


SEE 2 ! 1 , 


n 2 *, 2 

' a±n + ~2m <3in ' 


2 

(0) + ~~~ u j. n (0) 


+ w in (°) (5) 


in which w^ Ii (0) is the radial displacement coefficient of the 
sheet, of hay i at x^ a 0. 

Appropriate changes of the subscripts i in equations (3) 

■k° (5) permit the application of the equations to each hay of the 
structure ■ 

Sing deformations . - The radial displacement at any point cp 
of a synaaetrically loaded circular ring can he expressed as the 
Fourier expansion (see pp. 208 and 209 of reference 4) 


'CO 


cos nep 


It can he shown hy the method' of virtual work (pp. 209 and 210 
of reference 4) that for inextensional deformation the radial dis- 
placement coefficient ( i v iii) r j_ n g for a ring of radius E and 
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constant moment of inertia 1^ ■ that is loaded "by the shear flows 

in bays i and i-1 and by an arbitrary set of symmetrically applied 
external forces is (fig* 1) 




+ C 


in 


(6) 


In equation (6) the first expression on the right-hand side repre- 
sents the part of the radial displacement coefficient due to tho 
correction shears only, whereas the second expression represents 
the part of the displacement coefficient due to the external loading 
and the elementary shears. Values of C in are given later for 

particular loadings. (See equations (23).) 

Continuity relationships .- The following expressions can be 
obtained from continuity considerations of the rings and sheet of 
bays i-1, i,,and i+1 (fig. l) : 


( 7 ) 

t ’in( L i) - o i + l,n (0 ! ' < 8 > 

“l-l,n( L l-l) " W°> (9) 

u in( L i) ■ u i + l,n (0) <“> 
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*i-l,n(°) = ( v i-l,n) rins 
w i-l,n( L i-l) = ( w m) ring 


> 


( 11 ) 


^ = Wring 
"inW * Wiring 


> 


( 12 ) 


w i+l,n(°> = ( w i+l,n) rlng 
w i + l,n( L i+l) " ( V i+ 2 ,n) ring 


> (13) 


Equations (7) to (10) are conditions of continuity of a and u 
across the "boundaries "between bays i -1 and i and between bays i 
and i+1. Equations (11) to (13) state that the radial deformations 
of the rings bounding bays i- 1 , i ; and i +1 are eq.ua 1 to the sheet 
deformations of these bays at the rings. Implicit in equations (11) 
to ( 13 ) is a statement of the continuity of v of the cylinder 
across the boundaries between bays. 

Pft^nyysnc a formula Substitution of the expressions for 

u in( z i)> <S) rtog (eluatioDB (3) 

to ( 6 ). respectively) in the continuity relationships (equations ( 7 ) 
to ( 13 )) yields the following seven simultaneous equations in 
•which n = 2 , 3 , h, . . . 
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nL.« .t 




nL.» 

- 0 


pj . j 

^i+l ,iS 0 ) ~ a in^ 0 ^ '*"^7^ 
“ln (0) ‘ "1-X,a (0) + 


■H+l, 


n' 


t (0) - u^JO) + -S-h 

1 £SRt ^ 


4r a la " = 0 


R 4 8 in ' a i-l,n . „ *%-l n 3 ll l-l 3 

' ' 11 + C j es “ 1Tr a 4 m n „ *"* "-■ - r ■-- ‘i 

EI i n(n 2 - l) 2 Gt i-1 ’ ^’l-l 

2t 2 2t 

11 Li-i . , « %-x , v 

+ 

R 4 a i-l,n " a i-2,n „ 

+ Cj.x n 

^ n(n 2 - l)" 


■®i 


JL 

si. 


> i 

R 4 a l+l,n ” a in nil . 

n a a - 


‘( n£ - x T 

+ _^ r 

2ER 


•> ^ 

nLj_ , 


n 2 L± 2 , , n 2 L, , rA a in - a,- 

i~a (o) + iu, (0 + 5— J£ 3 

2ER ^ R in v/ ei 

A n^ - 

R 4 a i+2,n " a i+l,n . „ nL i+l 

ttt , + c i+2>n =* a i+l,r 

SI i+ s n(n 2 - l)- “l+l 


2 T 2 2 t 
n L4 + x n L. 

+ “S-CT-.! ^(0) + : 

2ER i+1 > n R 

R 4 a i+I,n " a ln 

EIi+1 n (n 2 - if 


f%kL 


■ U i+l ; n^ 0 ^ 


> 


1+1,11 


+ 


I+l,n 
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If the six quantities ff ± _ 1#a ( 0 ), ar^O), 

n (0) , u^ n (0), and 1^4.1 n (0) are eliminated from the 

seven expressions of equations (lU ) , the following recurrence 
formula relating the Foirier coefficients a of five successive 
bays is obtained: 




+ a 


i+1, n 


fe- .. , ha , *w - 

%+i x i+i\ x i+2 SAj.+ir 


U- 


- a.. 


i+2,n 


fi+S 


K°i-l, n • ( V 1 + ”2)°to + (®2 + ”3) 0 i + l,n. - ” 3 C i 4 - 2 , J-f- <«) 


in which 


». - _i Jbsi + 


1 L i-1 2 ^ Ll ^ +1 




= 


1 L i L i-i X 

2 \% + i + x ± + i + y 


U 3 


-4ft? 

J i+1 \ 1 / 
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» t . + 5 a=i + £1=1 + i 


4 Lj , 2 \Ll + l ■ L i+i ^ 


A i = 


Vi 3 


B i - 


St 

Gt i L i 2 


(n 2 - l) 


If the cylinder is of uniform construction, equation (15) can 
be considerably simplified and reduces to 


*1-2, n + 2 ?n a i-l,n + 2 l 3 n a in + 2r n a i+l,n + a i+2,n 


= ( C i-l,n " 3C ln + 3 °i+l,n " C i+2,n)“ 


in which 


„ n 2 - 6B 

7 = - 2 + 

n 12A7 


P„ = 3 + 


n 2 + 3B 
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The recurrence formulas ( 15 ) and ( 16 ) relats "the nth shear- 
flow coefficient pf "bay i with the corresponding coefficients of 
the two hays on each side of b^y i. .One equation similar to equa- 
tion (15) or equation (l6)‘ can ' consequently be written for each 
bay of a cylinder, provided that at least two bays exist on each 
side of this bay* For anti symmetrical loading, equations (15) 
and ( 16 ) can be applied if the- Fourier coefficients a are replaced 
by the coefficients b. 


Boundary Equations 

r 

Since the recurrence formula applies only to a bay having two 
bays on each side, incomplete or boundary equations must be found 
for each of the two bays at each boundary. Boundary equations, 
consequently, are presented for bays m and m-1 for a cylinder fixed 
at the right of bay m and for bays 0 and 1 for a cylinder free at 
the left end of bay 0 • (See fig. 3*) By suitable combinations of 
.the boundary equations and by proper manipulation of the subscripts, 
these equations can be used for the analysis of cylinders fixed at 
both ends, unrestrairiod at both ends, or unrestrained at one end 
and fixed at the other end. ‘ • 

Procedure for deriving boundary ' equations . - The general' recur- 
rence formula was derived by combining the equations for , 

l1 ln( x i ) > w ln.( x l) > and ( w in) ± (equations (3) to (6)) with the 

general. .‘continuity conditions (equations (7) to ( 13 )) and then 
eliminating all Fourier coefficients except the a's. In the deriva- 
tion of the boundary equations,' the .defining equations (3) to (6) 
are combined in a similar fashion -with (l) all of the continuity 
conditions (equations ( 7 ) to ( 13 )) that do not include quantities 
in nonexistant bays or rings and (2) the boundary conditions. 

Thus, for clamped edges (see fig. 3 ) the boundary equation for 
bay m is obtained by combining equations ( 3 ) to (6) with the con- 
tinuity conditions 


<r m-l,n( L m-l) a mn^ 0 ^ 

l Wl,n( :i i2i-l) = u mn(°) 

v m-l,n (o ) = (%-l,n) ring 
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Vi,d(Vi) = { w ““) rln ^ 

and the boundary conditions 

D 

' arid' -than eliminating all the Courier coefficients' except the a 'a. 

■ ■' Boundary equations for fixed end .-' If the foregoing- procedure 
is followed, the boundary actuation for bay m (fig. 3) is found to 
be. ......... 



in which 


# ■ 
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and similarly for lay m-1 



+ (>*5 + » 6 ) C J 4 - 


(18) 
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•where 




1 Am-2 
W * V '”" 1 




-JU&i + + ^ 

Vi 2 V ^ ^ / 




1 An-a 
" V \Vi + 




1 An-1 

Vi 2 U 


• L m~2 

*ta 



For cylinders of uniform construction the fixed-end boundary 
equations (17) and (l8) for bays m and m-JL, respectively, reduce to 


Ml 


'bl o „ + (Zy - lW , + (2[i - 27 - 6la 

m- 2 ,n \ n } m-l,n V, n. 'n. / : 


^.o + 2 VW,n + 2 &nVl,n + ( 2r n + ^'Wi 


> (19) 


- («n-2,«i ’ 3 c m-l, n + ^Vfi- 

'Irny 
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For anti symmetrical loading the Fourier coefficients a in 
equations (17) to (19) are replaced "by the corresporuiing coeffi- 
cients b . 

In order to apply equations (17) to (19) to the left end of a 
cylinder, the signs of the shear-flow coefficients must he changed 
and the subscripts of the various terms altered* If the cylinder 
of figure 3 is fixed at the left of hay 0, subscripts m, m-1, . • • 
are replaced hy 0,1, . . respectively, for those terms pertaining 
to the sheet of the hays and hy 1,2, . . ■ respectively, for those 
terms pertaining to the rings. 

Boundary equations for unrestrained end .- The boundary equa- 
tions for the unrestrained end of the cylinder shown in figure 3 
aro also found hy following the procedure outlined. The boundary 
condition at the free edge is 


a OrS 0) * 0 


The boundary equation for hay 0 is found to he 


"On 


\l f 

-ill + 

% V I; 


% 


6Bq - it 
6A 0 7 


Xo 


X 3 n 

=w 





re 


fx 


iTCa “ ( X 1 + * 2 )° 


In + x 2°2n 


1 E 




(20) 
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In ■which 


Similarly 
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L 0 L 1 



the "boundary equation for "bay 1 is 




( 21 ) 
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•where 


H 


i 

L l L 2 


, , iff* + y 

5 x 0 2 \ 1 i 


X, . — ( — + 1 

6 b. \ L 1 


Por cylinders of uniform construction the unrestrained-end 
boundary equations (20) and (21) for hays 0 and 1 , respectively, 
are 


( 27 n + 2P n + ^) a 0n + ( 27 n + X ) a ln + °2n = ’ ( 


n \.°0n ~ 2C ln + C ^nj h 


»): 


El 


ITny I 


a Von + £| 3 n a ln + ^n a 2 n + a 3 n ' (°Qn * + 30 Sn ' ° 3 n) 


El 

B^ny 


Por. antiBymnetrical loading the coefficients a are replaced 
"by b in equations (20) to (22) , In order to apply equations (20) 
to (22) to the right end of a cylinder, the signs of the shear- 
flow coefficients must be changed and the subscripts of the various 
terms suitably altered. 

Special boundary equations .- The boundary equations developed 
are suitable for cylinders having four or more bays. Por the 
special case of the center bay of a three bay cylinder, the boundary 
equation., "which depends upon the conditions at both boundaries, can 
also be found by means of the general procedure previously outlined. 
The boundary equations for cylinders of one or two bays can be 
similarly derived. 
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Application of Recurrence Formula and Boundary Equations 

Specific loadings .- As mentioned previously, the stress 
analysis of a reinforced cylinder arbitrarily loaded in the planes 
of the rings can he carried out conveniently if the stresses caused 
by the symmetric and antisymmetric components of the external 
forces are suitably combined. Further simplification of -the 
analysis is obtained if the loadings are resolved into concentrated 
radial forces, concentrated tangential forces, and concentrated 
bending moments. For each ring loaded at 9=0° (see fig. the 
load function C h ri , obtained in the derivation of equation (6) 

for ^ , for a concentrated radial force, a concentrated 

tangential force, and a concentrated bending moment are, respec- 
tively, 



P l n R^ny 

iB Eli 



^i R 4 ny 
j © El. 




> (23) 


where P, T, and M c are the symmetrical radial load, the anti- 

symmetrical tangential load, and the anti symmetrical bending-moment 
load, respectively, acting, on any ring i at cp = 0°. 

Simultaneous equations ■ - A typical set of equations applicable 
to a cantilevered uniform cylinder with six bays (m = 5 in- fig- 3) 
is presented in table 1. The first two and last two rows wore 
obtained from, the unrestrained- end and fixed-end boundary relation- 
ships of equations (22) and (19), respectively, end tho inter- 
mediate-rows were obtained from the recurrence formula of equa- 
tion (16) • For a nonuniform' cylinder these expressions are 
replaced by those of equations (20), (ul) , (lG), (1-7); and (15) . 

Xt is to be noted that tho coefficients of the unhncvn a 's 
and b’s are independent of (load term of equations (23)); 
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consequently, numerical solution of the equations (reference 5) 
for various loadings is greatly facilitated. A sot of simultaneous 
linear equations simi3ar to that of table 1 must he solved for 
each n-value chosen. The number of n-values, required depends upon 
the desired accuracy. The Fourier coefficients obtained for a 
given load P, 1, or M at cp = 0° can be used to determine 

the coefficionts for similar loads at any other value of cp since 
the z-axLs (fig. 2) can be chosen to coincide -with any radius. 

Stresses and loads in cylinder .- After the coefficients a 
and b are computed, substitution in the formulas (Al) to (A4) 

' presented, in- appendix A enables the stress analyst to compute the 
shear flow in the actual sheet, the direct stress in the fictitious 
sheet, and the moments, shears, and axial forces in the rings. The 
stresses due to loads acting at several rings and at various values 
of cp can be superimposed to give the stressos caused by these 
loads acting simultanoously . 


APPROXIMATE METHOD OF ANALYSIS BY SOLUTION 
OF FINITE DIFFERENCE EQUATION 
Difference -Equation Solution for Infinitely Long Cylinders 

Equation (15) referred to previously as a general recurrence 
formula is also a fourth-order finite difforence equation "with 
variable coefficients. Since the variable coefficients prohibit 
the solution of this equation in closed form, Only the solution of 
the equation that pertains to a uniform cylinder Is discussed 
herein. A general procedure for solving the fourth-order finite 
difference equation with constant coefficients (see equation (l6)) 
is presented in reference 6. When the right-hand side of equa- 
tion ( 16 ) Is sot equal to zero, the following homogeneous equation 
is obtained: 


a i-2,n + 2 ?n a I-l,n + 2 Pn a In + 2 VH+l,n + a i+£,n = 0 


(24) 


From reference 
consists of the 



6, the general solution of this homogeneous equation 
following six Independent solutions: for 

> 1 and < 0 
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for 


for 


for 


for 


for 


a in * ^(“Ln 003 k V °Sn aln k *n) 
lc 

+ e n (% cos k % n + % n Bln k t n ) ( £ 5a) 

D > 1 and 7 n > 0 

«»«!» 1c " 

a in = (_l)ke ^ (“in COS k X n + °2n sln k %n) 

+ cos k% n + a kn sin k% n ) ( ? 5h) 


D n < 1 and 7 n < 


-ii/ k 


a in * e n cosh fc P n + “211 BjLnh ^n) 

\j? k 

+ e n ^oc-jh cosh kp n + cc^ sinh kp^ 
D n < 1 and 7 n > 0 


a. 


in 


+ sli* ip n ) 


D n = 1 ®nd 7 n < 


-ll] 




Y J y. 

in = s (“m + < W C ) * « n («3n + %7) 


I> n = 1 and 7 n > 0 


(25c) 


= (-I) 11 © cos k kp n 

+ (-l) k e Vn (ctgn cosh kp n + o^ n sinh kp n ) ( 25 d) 


(250) 


a 


in 


= (-1) V'^'c^ + ^n k ) + (•' 1 ) k ^ nk (a 3n + o^k) (£5f) 
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k = i = 0,1,2 . . . 


and <^ n } and a^ n are arbitrary constants. 

The analysis of a uniform cylinder that extends longitudinally 
to infinity in both directions from a loadod ring is readily carried 
out with the aid of equations (l 6), ( 2 . 1 ), and (25), If the loaded 
ring is considered to be a boundary between the two halves of the 
beam and if no load other than that at the boundary is ass ume d to 
act, the difference equation (16) with tho right-hand side set 
equal to zero applies equally well to both parts of the cylinder 
(see fig. 5); consequently, only one-hulf of the cylinder need be 
considered in the analysis. Since the difference aquation 
applicable is the homogeneous equation (21), equations (25) 
together with the appropriate boundary conditions are solutions of 
the present problem. 

The distortions caused by the concentrated load have no effect 
On the stress distribution in the cylinder at k = 00 ; therefore, 

V = °* Ihe first term on the right-hand side of each of equa- 
tions (25) satisfies this condition^ however, the .second term does 
not satisfy this requirement and, hence, must vanish. The solutions 
then that are compatible with the boundary condition at infini ty 
are from equations (25)? for D n > 1 .and 7 n < 0 

a in = 6 ^(“ln cos k 'X h + «2n sln 


(26a) 
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for D n > 1 and 7 n > 0 

a in = ^ 1 ^ e " Vnk ( a ln cos *\ + “211 sin k ^) ( 26b ) 

for D n < 1 and 7 n < 0 

a in = e "^ nk ( ct ln cosh kp n + ^n einh kp n) ( 2 6c) 

for D n < 1 and 7 n > 0 

, -\j/ k . 

a ln = e n ( ct in cosh kp n + “Sn Billh ^n) ( 2 6d) 

for D n = 1 and 7 n < 0 

a ln “ 0 ^(* 2 X 1 + °2n k ) ( 26e ) 

for D_ = 1 and y_ > 0 
n n 

a ln = (_l)ke ^Kn + ^n 1 ") < 26f > 

From the conditions of symmetry about the loaded ring, 
modification of equation (l6) leads to the determination of two 
boundary equations applicable to the present problem. If the 
load function at the loaded ring is designated CQ n (see equa- 
tions (23)) and equation (16) is written for bay 0, the first 
boundary equation is ; 

( S|3 n ‘ ^nKn + (^n ' + «2n = -3 C 0n 

E^n 7 


(27) 
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since a~ = -a n and a-,„ = -a 0 . If equation (l6) is written 
On -In -i*n -^n 

for bay 1, the second boundary equation is seen to be 


(r. 


( 28 ) 


since a» = -a n - 
On -In 


The boundary equations ( 27 ) and (28) permit the determination 
of the arbitrary constants cc-^ and -^ or a Si ven value of n, 

substitution of the appropriate value for a^ n from equations (26) 

into equations ( 27 ) and ( 28 ) yields a sot of two simultaneous equa- 
tions; for example, if I> n > 1 and 7 n < 0 


<X ln 


-\|/ -2 \{/ 

( 2p n " 27 n) + ( 27 n ' 1 ) e * COS \ + 6 “ 008 


°2n 


-i|/ -2i ! 

(2/ n - l)e n sin % n + e n sin 2% n 


= -3C, 


El 


0n T,It. 

BV 


a ln 


_\i/ -3*4/ 

|(£7 n - +2P n e 31 cos 2/^e 11 cos 2^ + e 11 cos 3^ 



sin % +2 ye 
A 'n n 




sin 2% 


£ V 0 





El 

En7 


The constants a 1r) and ctg^ are obtained from the solution of 

these equations. To each value of n there corresponds one value 
each for a 1y , and ctrp ^ . Since D n and 7 n are functions of n 1 

well as the elastic properties of the cylinder, for a particular 
cylinder more than one of equations (26) may be required, for the 
determination of all the values of and cv^. With the values 

of these constants determined for each value of n, corresponding 
rallies of a in for each bay are obtained from equations (26) . 
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As In the application of the recurrence formula, a’s corre- 
sponding to several harmonics, that is, n varying from 2 to the 
value that yields the desired accuracy, must be found# For anti- 
symmetrical loading, a is replaced by b in equations (24) 
to (29) • The values of the coefficients a and b obtained are 
substituted in equations (Al) , (A2) , and (A4) for tho desired load 
values. Since the expression for the direct stress in the sheet 
now involves an infinite summation along the cylinder of the shear- 
flow coefficients, simplified formulas for the direct stresses at 
any ring k are presented in appendix B. 

If equations (27) and (28) are replaced by the unrestrained-end 
boundary equations (22), with all values of C_, n except sot 

equal to zero, a tip loaded cylinder extending to infinity in one 
direction can be analyzed with a procedure similar to that developed 
herein. 


Application to Finite Cylinders 

■Whereas a concentrated load causes distortion in the region 
in the immediate vicinity of the load, for most practical purposes 
the part of tho cylinder a few bays away from the load can be 
assumed undisturbed. Consequently, if the load is located a 
sufficient distance from external restraints, the distortions of 
the cylinder in the region of the load, are independent of those 
restraints. If then a uniform cyddnder of finite length Is to bo 
analyzed and. this cylinder Is loaded in a manner auch that the 
load is not in the proximity of an external restraint, the ele- 
mentary stresses and loads are found as usual by considering the 
cylinder to be finite, whereas the corrections may be found by 
use of the difference-equation method by considering the cylinder 
to be infinite. Because the effoct of the concentrated load 
dissipates quite rapidly, values of a and b are usually of 
interest only for those bays in the vicinity of the load. The 
desired forces and moments in this region can then be determined 
as before from the equations given in appendixes A and B. 


Adequacy of Difference -Equation Solution 

Although the solution in closed form of the problem of a 
uniform reinforced circular cylinder is exact only for infinitely 
long cylinders symmetrical about a loaded ring, comparisons of the 
finite-difference -equation solution, tho re curronco -formula solution, 
the standard solution (reference 7)/ and experimental data for 
cylinder 2 of reference 3 were made for a cylinder fixed at one 
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end, unrestrained at' the otlie!r, and' having oi3^ ; *f*b^ir Tjayb^ .The. 
cylinder-.- s®s loaded with a- concentrated radial; force It; a ring 
located '.two* tays:‘ firm each ' end . ' '■ ( See '■ fig. ;■&*)■* • ‘in' figured; J‘ aijd 8 
curves '■ are ".given • foir-' bending ■ moment’s in - the -loade’^I ring and^ad jacent 

« ■ . *1 *1 -a j.i_ _ mn 1 ’ ' r 2-. - — JM ■& — >k-u'J« 1 4- a 



few bays:; an extreme case is represented, that 'is unlikely to b.e' 
met. .in practice . .'Ihe' more' bays a cylinder IiaS' the’ inbre clo'sely 
it approximates an infinite cylinder for which ilie finite -diff epehce- 
equatlo.n. solution is . exact; ■‘cO'ftsequeritiyy the; gc'od .Agreement shown 
in f igitres _T . and 8 among the; 'finite -d’ifferehce -equation solution, ■ 
the recurrence-formula solution and" experimental ' cLa t*J' indicates, 
that the; simplified solution is-.'qutte adequate.' * ' ' ", 


Advantages of Difference -Eq.ua tion Solution 


Since airplane fuselages ' approximating' circular' cylinders are v 
composed of a relatively- -large numb er : of -bays ;.-f or' moist practical 
cases, the simplified solution- should be a good approximation to 
that obtained by the use of the recurrence formula. As mentioned . 
previously, when the recurrence formula is applied, sets of 
simultaneous equations containing as many unknowns as there are 
bays in the structure must be solved for each n-value required. 

For structures having many bays the amount of computations involved 
may be prohibitive; however, no such computations are involved when 
use is made of the infinite -cylinder solution. In addition, this 
solution is adaptable to the construction of design charts similar 
to Wise 's charts of reference 7 • The analysis of any long uniform 
cylinder is dependent only on the values of the structural 
parameters A and A/B. For various representative values of 
these parameters, charts can be constructed from which the analyst 
can determine desired stress coefficients. For extreme cases, such 
as a cylinder loaded only one to two bays away from a restraint, 
the recurrence-formula method is recommended for accurate solutions.' 


COITCLUDING REMARKS 


The recurrence formula developed in the present paper facili- 
tates the stress analysis of circular cylinders loaded in the planes 
of the reinforcing rings. The cylinders can be composed of bays 
of different cross sections and lengths and can be supported by 
rings having different moments of inertia. The boundary equations 
presented are applicable to cylinders fixed at both ends, unrestrained 
at both ends, or unrestrained at one end and fixed at the other end. 
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Tor the. analysis of cylinders composed of Relatively: few 
bays, it'ie recommended that 'the recurrence formula ■■'be used to, ' 
obtain sets .'of simultaneous linear algebraic .equations.' The-., 
solutions of these equations lead'-tccian accurate- determination- of 
the stresses- in the rings . and sheet of -the cylinders. The analysis 
of cylinders : composed of many bays, as are semimonocoque fuselages, 
can more conveniently be accomplished by the solution of the 
recurrence 'formula.'asa finite -difference equation. .-Although the 
stresses • obtained with this .ablution. are approximations to -the 
more: accurate . stresses found with the simultaneous .equations, for 
long. cylinders the, computations- involved 'are considerably shorter. 

In addition,' -.since for the three basic loads the. stresses .determined 
by this method are dependent- only upon the structural .'parameters 
of the cylinder, charts facilitating the rapid determination of 
the stresses in reinforced cylinders can be readily constructed. 


Langley Memorial Aeronautical .liSbaratory •’ 

.. National Advisory Committee for Aeroiiautics ' .. 

• ' Langley Field, Va.y November 12, 1946 i 
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APPENDIX A 


FORMULAS FOB LOADS AND STRESSES IN CYLINDERS 


After the coefficients a and b are computed, the shear 
flows in the actual sheet, the direct stresses in the fictitious 
sheet, and the bending moments, shears, and axial forces in the 
rings can be found -with the aid of the equations • given in the 
appendix of reference 3. For the sake of completeness these eq.ua 
tions, -with some additions, are presented herein. 


Shear Flow 

The total shear flow q^(<p) in any bay i for any ring 
loading on a cylinder aan be expressed as 




% 


CD 

* 1 . 

n=2 


in 


sin ncp + 


00 

I* 

n=2 


cob ncp 


(Al) 


in which q^ represents the elementary shear flow calculated on 

the basis of rigid rings. For a cantilevered cylinder, q^ is 

zero for those bays located between the tip and a loaded ring. For 
those bays between a loaded ring and the root, tho values of 

for a radial load P, a tangential load T, and a concentrated 
ring bending moment M Q , each applied to ring i at 9 = 0°, are 

given in table 2. Positive forces and bending moments are indi- 
cated in figure 4, If more than one ring is loaded or if the 
cylinder is not of cantilever construction, P^, 11, and M c ^ 

are replaced by the resultant radial, tangential, and moment load, 
respectively acting on a cross section of bay.i. 


Direct Stress in Skin 

t 

For a cantilevered- cylinder such as that shown in figure 3, 
if the longitudinal skin stress at ring 0 is assumed to be zero, 
the direct stress at ring i is (see equations (lb) and (3)) 
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0.(0, cp) = 0-n — ) la. + a,„ + . 

i Y R R L— \0n t’ 111 t % 


L i-1 


i-l,n t'. 


n cos rxp 


?On^" + t ln^T- + • * ( + Vl ; np — )nsin n<p (A2) 

V r 0 ■ 1 ' t i-l/ 


in which 0 E is the stress given by the simple engineering theory 

of bending. Since the shear stress is constant in the longitudinal 
direction -within a bay, a varies linearly between rings. 

If the cylinder is rigidly fixed at ring 0 .as well as at 


ring m+1, the initial boundary stress y Cq^u; cob ncp i^ror 

n=2 

symmetrical loads) must be added to the direct stress obtained with 
equation (A2) . The value of the Fourier coefficient ^on^ is 

determined for a cylinder having at least three bays from the con- 
tinuity condition 


0 n (0) cos ntp (for 


w On( L o) = ~( w ln) . 

/ v 'ring 


and the boundary conditions 


UQ n (°) = 0 


w tVi(0) = 0 


together with the defining equations (5) and (6). The relationship 
obtained is 


2Lo A ' 0 7 f <SB 0 - n 2 \ El! 

"°n <0) = - “onf 1 + 7 V + nr°ln (A3) 

0 \ 0 7 J ITny 
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in which 

B S t' 0 

Vo 3 

For anti symmetrical loading, a and a are replaced 

On In 

"by tg n and b^, respectively. 


Bending Moments and Forces in Ring 

The bending moments, shear forces, and axial forces in the 
reinforcing rings of a cylinder arbitrarily supported at its ends 
are, respectively. 



“i 


+ R 2 

n=£ 11 ( n2 " l) 


cos ncp 




in which Mg, and Hg are the bending moment, shear force, 

and axial force in the rings, respectively, determined on the basis 
of elementary shear flow in the skin. Positive values of the bending 
moments and loads in a cylinder are indicated in figure k. Form ula s 
for Mg, Y R3 - and -.Hg corresponding to a radial load P, a tan- 
gential load T, and a concentrated ring bending moment M c , each 

applied to a ring i at <p = 0°, are given in table 2. For rings 
not loaded externally, only the series expression in equations (Ak) 
are required. 
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APPENDIX B 


DIRECT STRESSES IN INFINITELY LONG CYLINDERS 


Eor the determination of the direct stress in the sk in of a 
uniform infinitely long cylinder, equation (A 2 ) can he replaced 
hy 


<^(0, cp) 


^R 


03 



n ='2 


k 

n cos ncp 

i— CO 


(Bl) 


or 


ojo, cp) 





8 in 




n cos ncp 


(Bl=) 


in which only the coefficients a^ n are considered. In these 
equations, P-fX 0 ; <?) is the direct stress at ring i=k. Corre- 
sponding to the six values of from equations (26), six 

solutions for 0^(0, 9) can "be determined ay summation along the 

cylinder. As an illustration of the procedure involved, equa- 
tion (26a) 1 b used herein for the value of a, . Consequently, 

equation (B 2 ) becomes 


cr^O, cp) = 0 B + 


Rt’ 


V" 

ZL 

n=2 


V 'V 


|Z_ 

iso 


Kn 


00s ^ X n ' I ’ a £ n sin. k 


'X.,) 


■ Z 8 “ (“in 005 Sin =:%:,) 

i=0 


n cos ncp (B3) 



The emanations from 1 = 0 to i = » and i = 0 to i = k ~ 1 are readily accomplished with 
the aid of formulas 6.830, 6.833, and 3.0I, numbers 12 and 13 , of reference 8. The resulting 
formula for the direct stress at any rin^ k is for >1, 7 B < 0, and i ■= k = 0,1,2, . . . 


00 


, . i V 

’’ )=0 e + 5^ .4. e — 

s ■ • ■ , n=2 


E* 


cos k^-cos (k-1) 


%n +a 2n *** 


sin kX^-sin (k-1) 




cosh - cos X n 


n cos ncp 

I 

m 


Vith a procedure analogous to that used for the dete rmina tion of this equation, the following 
solutions are obtained for the direct stresses corresponding to the remaining five values of a in : 

for D n > 1 and 7^ > 0 


00 


yi 
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{^(O.cp) 

for' ft 
n 

^(o.q?) 

for D 
n 

for D n 

(Tl(0,{p) 

for D 

n 

<^(0,9) 

If 

to (B5) 


r.v 


i 

9 


tr, 


, k L V -* n * °3n [ e C0B k V «»0^ 1 )^ n J+ j* ^ *V eln(t-l)^| 

: Op+ (-1) - — / & n — — - n cos rap (B5a) 

iEt’— cosh^ cos 


CO 

o\ 


rt=^ 


< 1 and 7 n < 0 


oo 


t 


_ j _. y 1 

x, je n cosh kp n - coah(k-l)p n J + a £n je ra slnli kp n - einh(k-l)p. 


1 Cfp + / © 

K URt 

n=£ 


d +0 4? 




J 


cosh ik - cosh p„ 
T n n 


— n cos’ rap (B5h) 


< 1 and 7 n > 0 




00 ^ k a j J1 e 1 ^ 1 cosh kp^+coshCk-lJp^ +^ 2 n e^sinh kp^+sir 


e 




n cos rap (B5c) 


rt=2 


cosh ij/ + cosh 


= 1 and 7 n < 0 


, L \ ' -%* 

v — ■ ■ -| — / 6 ■■ ■■ ■ 


( e^ 1 - l) + a^L'k - (k - 1)1 


= da -v / e 

* £Rt • l— 


H=E 


cosh ilr - 1 
■n 


el cob rap 


(B5d) 


= 1 and 7 n >0 


/ _,k L V A k ^ 

“ «*D + (- 1 ) — / 6 

cr>4- i / 


(e V:i + l) + a^Jke^ 11 + (k - 1)] 


2 Rt 1 L- 

Ttzd 


n cob rap 


cosh y n + 1 


(B5e) 


the coefficients b^ are considered,, cos rap is replaced by -sin 09 in equations (Bl) 
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Figure 3 ~ Side view of cantilevered cylinder. 
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Figure 4.~ Sign convention used in analysis 
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Figure 5. - Loaded part of infinitely long cylinder. 
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Figure 6 “Side view of cylinder 2 analyzed in reference 3. 
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Fig. 7 
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Figure 7.- Comparison between calculated and experimental ring-bending 

moments for cantilevered cylinder. 












